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ABSTRACT 

We apply standard disk formation theory with adiabatic contraction within cuspy halo models 
predicted by the standard ACDM cosmology. The resulting models are confronted with the broad 
range of observational data available for the Milky Way and M31 galaxies. We find that there 
is a narrow range of parameters which can satisfy the observational constraints, but within this 
range, the models score remarkably well. Our favored models have virial masses 10 12 M a and 
1.6 x 10 12 M for the Galaxy and for M31 respectively, average spin parameters A w 0.03 — 0.05, 
and concentrations C vlr = 10—17, typical for halos of this mass in the standard ACDM cosmology. 
The models require neither dark matter modifications nor flat cores to fit the observational data. 
We explore two types of models, with and without the exchange of angular momentum between 
the dark matter and the baryons. The models without exchange give reasonable rotation curves, 
fulfill constraints in the solar neighborhood, and satisfy constraints at larger radii, but they 
may be problematic for fast rotating central bars. We explore models in which the baryons 
experience additional contraction due to loss of angular momentum to the surrounding dark 
matter (perhaps via a bar- like mode). These models produce similar global properties, but the 
dark matter is only a fourth of the total mass in the central 3 kpc region, allowing a fast rotating 
bar to persist. According to preliminary calculations, our model galaxies probably have sufficient 
baryonic mass in the central ~ 3.5 kpc to reproduce recent observational values of the optical 
depth to microlensing events towards the Galactic center. Our dynamical models unequivocally 
require that about half of all the gas inside the virial radius must not be in the disk or in the bulge, 
a result that is obtained naturally in standard semi-analytic models. Assuming that the Milky 
Way is "typical", we investigate whether the range of virial masses allowed by our dynamical 
models is compatible with constraints from the galaxy luminosity function. We find that if the 
Milky Way has a luminosity Mk — —24.0, then these constraints are satisfied, but if it is more 
luminous (as expected if it lies on the Tully-Fisher relation), then the predicted space density is 
larger than the observed space density of galaxies of the corresponding luminosity by a factor of 
1.5-2. We conclude that observed rotation curves and dynamical properties of "normal" spiral 
galaxies appear to be consistent with standard ACDM. 

Subject headings: cosmology:theory - galaxy structure 

1. Introduction sical problem which has seen many past itera- 

tions (Einasto 1972; Innanen 1973; Schmidt 1975; 
Modeling the mass distribution within the Einasto 1979). Because of our unique ability to ob- 

Milky Way (MW) and M31 galaxies is a clas- 
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serve these systems in exquisite detail, this prob- 
lem has continued to be a fertile testing ground for 
theories of galactic structure and galaxy formation 
(Kent 1989; Dehnen & Binncy 1998; Wilkinson & 
Evans 1999; Evans & Wilkinson 2000; Oiling & 
Merrifield 2001). In the absence of a coherent 
halo formation theory, early models were plagued 
by the arbitrary structure parameters of the dark 
halos. Recent developments in the theory of dark 
matter halo formation make this sort of modeling 
much more meaningful because the halo parame- 
ters (such as total mass, peak rotation speed and 
concentration parameter) may no longer be con- 
sidered free parameters but are correlated, and 
these correlations are specified within the con- 
text of a given Cold Dark Matter cosmology (e.g. 
Bullock et al. 2001, and references therein). 

While very successful on large scales, the stan- 
dard cosmological theory is experiencing mount- 
ing difficulties on the scales of galaxies and dwarf 
galaxies. The so-called "sub-structure" problem 
(the predicted over-abundance of sub-halos rela- 
tive to observed dwarf satellites of the MW and 
the M31 galaxies (Klypin ct al. 1999; Moore et al. 
1999)) may have been successfully resolved by in- 
voking the presence of a photoionizing background 
that "squelches" star-formation in small halos 
(Bullock, Kravtsov & Weinberg 2000; Somerville 
2001; Benson et al. 2001). However, the so- 
called "cusp" and "concentration" problems re- 
main thorny. The most acute manifestations of 
these problems are encountered for objects which 
are dark matter dominated — namely, the cen- 
tral parts of low surface brightness galaxies and 
dwarf galaxies (Moore 1994; Flores & Primack 
1994; Moore et al. 1999; van den Bosch & Swa- 
ters 2000). For example, the high central density 
of ACDM dark matter halos appears to be in- 
consistent with the observed rotation curves of 
dwarf and low surface brightness galaxies (Blais- 
Ouellette et al. 1999; Cote, Carignan, & Freeman 
2000; Blais-Oucllctte, Amram, & Carignan 2001; 
de Blok et al. 2001; Alam, Bullock, & Weinberg 
2001). 

It is interesting to ask how ACDM models fare 
in their predictions for luminous, high surface 
brightness galaxies, such as the Milky Way and 
M31. These systems probe the dark matter prop- 
erties on mass scales of 10 12 M©, which are signifi- 
cantly larger than the 10 8 — 10 9 M scales tested by 



inner rotation curves of dwarf galaxies. The Milky 
Way is particularly interesting for constraining the 
central cusp of the dark matter because the inner 
3-10 kpc region of the Milky Way is well studied. 
Apart from the rotation curve and integrated light 
distribution, the Milky Way offers many more op- 
portunities to study very detailed data which are 
not generally available for extragalactic systems, 
e.g. individual stellar radial velocities and proper 
motions in the solar neighbourhood, and recently, 
Galactic bulge microlensing events. 

It is likely that the halo of our Galaxy con- 
sists mostly of weakly interacting elementary par- 
ticles (wimps) rather than black holes or failed 
stars (machos) , a necessary condition to apply the 
standard cosmological approach. The very low mi- 
crolensing counts by machos in the inner 50 kpc 
halo of the MW imply that elementary particles 
make up > (80 - 90)% of the halo (Alcock et al. 
2000; Lasserre et al. 2000; Zhao & Evans 2000). 
Therefore it should be valid to apply halo forma- 
tion models to our Galaxy, which predict an outer 
dark matter profile of the form p oc r~ 3 . In this 
regard, earlier halo models using cored isothermal 
halos with a p oc r~ 2 outer profile and a divergent 
mass (e.g. Oiling & Merrifield 2001), or very steep 
outer profiles p oc r~ 5 (e.g. Evans & Wilkinson 
2000) are unphysical and should be updated. 

Navarro & Steinmetz (2000) and Eke, Navarro, 
& Steinmetz (2001) were the first to use mass mod- 
eling of the Milky Way to constrain cosmological 
models. The dark mass inside the solar radius 
was compared with observational limits derived 
from previous mass modeling. This could give mis- 
leading results because of inconsistent modeling of 
the halo profiles. Also unlike in galaxies with low 
baryon content such as dwarf or LSB galaxies, the 
density of the inner halo can deviate significantly 
from the original Navarro et al. (NFW; 1996) pro- 
file in high-surface brightness galaxies because of 
strong compression during the baryonic collapse 
(Mo, Mao, & White 1998). Given the coupling of 
the inner and outer parts of the halo, it is neces- 
sary to repeat the mass modeling of the Milky Way 
and M31 including a self-consistent treatment of 
the baryons and dark matter component. 

The amount of dark matter inside the central 
~ 3 kpc of the Milky Way and M31 introduces 
another layer of problems. It is very likely that 
both galaxies have fast rotating bars. A rotat- 
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ing bar experiences dynamical friction induced by 
the dark matter. If the mass of the dark matter 
surrounding the bar is too large, the dynamical 
friction could slow down the bar rotation over a 
short timescale (Weinberg 1985; Debattista & Sell- 
wood 2000), in which case, we should not expect 
to see fast rotating bars. This argument some- 
times is presented as the most severe problem the 
ACDM cosmological models encounter on small 
scales (Sellwood 2000; Sellwood & Kosowski 2001; 
Evans 2001; Binney & Evans 2001). It is clear that 
some effect of this kind must exist, but a realistic 
estimate of its magnitude is still difficult for cur- 
rent analytic works (Weinberg 1985) and N— body 
simulations (Debattista & Sellwood 2000). 

The COBE map of the Galaxy and the high 
rate of microlensing events towards the Galactic 
bulge have been used as arguments that most of 
the inner Galaxy's mass is in a heavy stellar bar 
(Zhao 1996; Zhao & Mao 1996). This leaves very 
little room for a massive dark matter halo (Bin- 
ney, Bissantz, & Gerhard 2000). We address these 
issues in this paper, and will treat them in more 
detail in a paper in preparation (Zhao et al. 2001). 

Modeling the Milky Way and M31 turns out to 
be challenging because of their complex formation 
processes. The present-day dark matter halos can 
have density profiles very different from the origi- 
nal dark matter prediction (e.g. the NFW profile). 
They are squeezed and deformed by the gravita- 
tional force of the baryons, which collapse out of 
the NFW halo to form heavy disks and dense nu- 
clear bulges. The amount of squeezing also de- 
pends on the angular momentum transfer between 
the baryons and the dark matter halo. The dark 
matter cusp is ellusive to detection because the 
gravitational force at the center is dominated by 
the baryons. Apart from detailed data, the Milky 
Way provides a good laboratory for studying the 
transfer of angular momentum from the rapidly 
rotating baryons to the dark matter halo through 
dynamical friction. Such a transfer must have 
taken place in the MW and M31 since both have 
a prominent triaxial bulge/bar and spiral arms in 
their disks. An analytic formulation for the angu- 
lar momentum transfer is given in an accompany- 
ing letter (Paper III). 

Studying the MW halo can set general limits 
on dark matter halos only if the MW is typical 
among galaxies that form in a dark matter halo of 



a given mass. This "typicality" hypothesis can be 
tested to some degree with our parallel study of 
the nearby galaxy M31. It can also be tested by 
producing many realizations of a Milky Way-like 
galaxy using a Monte Carlo-based semi-analytic 
model of galaxy formation. The semi-analytic ap- 
proach allows us to make statistical predictions 
of the formation histories of ensembles of galaxies 
hosted by halos of a given mass, and to incorpo- 
rate the combined stochastic effects of distribu- 
tions of mass accretion and merging history, halo 
spin parameter, and gas accretion and star forma- 
tion history. This approach has the unique ability 
to allow us to estimate the chances of finding a 
certain combination of structural parameters that 
simultaneously produce the observed stellar popu- 
lations and gravitational potential in galaxies like 
the MW and M31. We pursue this question in a 
companion paper (Paper II). 

In summary, the goal of this series of papers 
is to test whether the current standard theory of 
galaxy formation, which seems to be very success- 
ful on large scales and very problematic on dwarf 
galaxy scales, can produce L* galaxies such as 
the MW and M31. In this paper, this is done 
by combining classical dynamical modeling of de- 
tailed data for the MW and M31 with current 
cosmologically-based models of dark halos and 
disk formation. 

2. Model Ingredients 
2.1. Mass components 

In traditional mass models of the MW one de- 
composes the baryons into several descriptive com- 
ponents: the nucleus, the bulge, the bar, the 
spheroid, the thin disk, the thick disk, and the 
cold interstellar medium in the disk. It is some- 
what problematic to determine where a compo- 
nent starts and ends, and this invariably adds to 
confusion of terminology. For example, the bar 
sits in the overlap region between the bulge and 
the disk, and most stars in this region rotate in the 
same sense with similar angular speed. So is the 
bar a distinct identity on top of an oblate bulge, 
is it merely the triaxial bulge, or is it the non- 
axisymmctric part of the disk inside corotation? 
The lack of a 3D map of the Milky Way makes it 
uncertain whether or not the disk is truncated in- 
side the corotation of the bar. These distinctions 
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are subjective at some level and sometimes unnec- 
essary when the important quantity is the total 
mass distribution of the baryons. 

We take a different approach. We divide up 
the mass of the Galaxy into only two components: 
a dark halo component and a summed-up bary- 
onic component. We describe this baryonic sum 
by a function M b (< r), which is the combined 
mass of baryons inside the radius r. It is un- 
derstood that the mass distribution is far from 
spherical. In fact it is highly flattened, and tri- 
axial inside ~ 3.5 kpc. We incorporate some mini- 
mal modeling of the dynamical effects of flattening 
and triaxiality, but without repeating earlier more 
rigorous 3-dimensional triaxial dynamical models 
(Zhao 1996; Hafner et al. 2000). 

Our models for the bulge/bar are motivated by 
the Zhao (1996) model of the galactic bulge/bar. 
That model used COBE DIRBE de-rcddcncd in- 
frared maps (Weiland et al. 1994) and a collec- 
tion of stellar kinematics data in the direction of 
the central parts of the galactic bulge. The fit to 
the data indicated that the bar has an elongated 
boxy triaxial shape. Zhao (1996) added a steep 
oblate nucleus to the best-fit so-called "G2" model 
of Dwck ct al. (1995). The nucleus and a massive 
central black hole are important for modeling the 
compressed density of the dark matter within 100 
pc. The total mass of the bulge, the bar, and the 
nucleus was (2.2 ± 0.2) x 10 10 M Q . No dark mat- 
ter was explicitly included in the model of Zhao 
(1996). A massive Miyamo-Nagai analytic disk 
potential was used to make the rotation curve flat 
within 8 kpc. Here, instead of the Miyamo-Nagai 
disk, we use a double-exponential disk as in Kent, 
Dame, & Fazio (1991). The combined density of 
nucleus, bulge/bar, and disk are modeled by the 
following three components: 
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where pi j0 , P2,o, and /9 3i o are characteristic densi- 
ties, determined by the corresponding total masses 
Mi, M 2 , and M 3 of the components pi, p 2 , and pz- 
The dimensionless radii s\, s 2 , and S3 are given by 
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The above aspect ratios in the components pi(sl), 
p 2 (s2) and ps(s3) are taken from Zhao (1996) to 
reproduce the shape of the COBE bar and the disk 
of Kent, Dame, & Fazio (1991). We take the ver- 
tical scale height zq = 400pc from Zhao (1996), 
and take the disk scale-length to be = 3000 
or 3500 pc, and the distance to the Galactic cen- 
ter to be Ro = 8000 or 8500 pc. These param- 
eter combinations produce acceptable fits to the 
light distribution of the Milky Way after proper 
normalization of the three components. Zhao 
(1996) normalizes the pi and p 2 components so 
that they have equal strength at about a distance 
z = 0.75zo on the minor axis. This fixes the 
ratio pifl ■ p2,o so that pi + p 2 joins smoothly 
from an observed power-law nucleus to the ob- 
served COBE bar. The ratio p 2 ,o : p3,o is fixed 
so that (Mi + Mi) : M 3 = 1:5, approximately 
matching the bulge-to-disk K-band luminosity ra- 
tio 1.1 x 10 10 L Q : 4.9 x 1O 1O L of Kent, Dame, 
& Fazio (1991). Kent et al. in fact introduced 
a linear tapering of their disk scale length from 
250 pc at the solar radius to 167 pc at the center 
for better matching of the observed surface bright- 
ness map. This makes their disk less massive than 
our double-exponential disk if both are normalized 
to the surface density and the volume density at 
the solar radius. The same is true of the triaxial 
models of Freudenreich (1998), where the axisym- 
metric disk is tailored into a triaxial hole inside 
3.5 kpc. 

We do not introduce this fine tailoring of the 
disk, instead favoring the mathematical simplicity 
of a pi bar and a p 3 disk, rather than a true bar 
or disk. The extra mass contained in our p 2 disk 
is at the price of reducing the mass of our p 3 bar 
given the same budget of mass or light within the 
central 3.5 kpc. 

We have deliberately chosen not to label the 
three components pi , p 2 and p 3 as nucleus, bar and 
disk given the uncertain truncations of these lat- 
ter components. Instead they should be treated as 
convenient functions to describe the overall distri- 
bution of baryons. There is no separate thick disk 



4 



component. When computing the rotation curve 
we approximate p 3 as a razor-thin disk, and pi 
and pi as spherical. The enclosed baryonic mass 
is computed by spherically averaging p b : 



M b {< r) = m BH + At: {p b )r z dr, (8) 
Jo 



where 



r dn 



(9) 



A more user-friendly analytic approximation of 
M b (r) is: 



M b (< r) = m BH + Mi [1 - cxp(-1.15xi) 



M 2 



+ 



erf (^)-^ 2eX P(-*) 
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where x\ = r/ri, x 2 = r/r 2 , x 3 = r/r^ and n 
and Mi are the scale-length and mass of the nu- 
cleus, r 2 and M 2 the bulge/bar, and r 3 and M 3 the 
disk. We use eq. (8)in our calculations. The total 
mass of the Galaxy inside a radius r is computed 
as 

M(r) = M h (r) + M dm (r), (11) 

where Md m (r) is the mass profile of the dark mat- 
ter. We include a central black hole of mass 
m BH = 2.7 x 10 6 M Q in the Milky Way models 
and m BH = 3.7 x 1O 7 M in the M31 models. 

In our Milky Way models the parameters of the 
bulge are kept close to those of the original Zhao 
(1996) bulge. This is done to preserve the agree- 
ment with the COBE data and with estimates of 
stellar radial velocities. Another reason is the mi- 
crolensing counts, which may be a problem for 
models of the Milky Way (Zhao & Mao 1996). The 
rate of microlensing events was estimated for the 
Zhao (1996) bulge and was found to be consistent 
with the constraints on the counts (Peale 1998). 

2.2. Dark matter halo properties 

A sparse set of observational data can generally 
be fit by a large range of models, and sometimes 
unphysical models give excellent fits. Dehnen & 
Binncy (1998) give plenty of examples of this kind. 
To obtain physically interesting models, it is im- 
portant to impose constraints based on a physi- 
cal theory of halo formation, and to select mod- 
els which are not only compatible with data but 



are also physically well- motivated. In this paper 
we use three results from halo formation theory: 
the shape of the mass density profile, the halo 
concentration-mass relationship, and the distribu- 
tion of halo spin parameter (angular momentum) . 

Throughout this paper, we assume that the 
Hubble constant is equal to H = 70km/s/Mpc, 
the Universe is flat and the contribution of matter 
(dark and baryonic) to the critical density is equal 
to fio = 0.3. We assume that initially the dark 
matter density profile is described by an NFW pro- 
file (Navarro et al. 1997, NFW): 

Ps 
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c = 

M vir = 
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where C and M v j r are the halo concentration and 
virial mass, and r v i r is the virial radius. In the 
above equations, the parameter p cr is the critical 
density of the Universe and 5th is the overden- 
sity of a collapsed object in the "top-hat" collapse 
model (5th ~ 340 for our cosmological model). 
Two independent parameters — C and M v ; r - 
completely define all relevant halo properties. 

For some halos the central slope may be even 
steeper (-1.5) than the -1 slope in eq. (12) (e.g., 
Moore et al. 1998; Ghigna et al. 2000; Klypin ct 
al. 2001). For radii larger than 0.5-1% of the virial 
radius (1-3 kpc for our galaxies) the difference be- 
tween this profile and eq. (12) is relatively small 
(Klypin et al. 2001). Baryonic infall has a ten- 
dency to reduce this difference even more, as we 
will discuss in section 4.2. Because the baryonic 
compression is very large and still quite uncertain 
and because of the relatively small differences in 
the density profiles, we will not consider profiles 
with a steeper -1.5 inner slope here. 

While we have introduced the halo concentra- 
tion C as an independent parameter, in practice it 
is found to be strongly correlated with the virial 
mass (Navarro et al. 1997; Bullock et al. 2001). 
For the ACDM model, Bullock et al. (2001) give 
C = 15 - 'S.SlogiM^/lO^h^MQ). Statisti- 
cal ensembles of halos at fixed mass show 40% 
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rms fluctuations of the concentration around this 
mean. For 10 12 M halos this gives a range of 
C = 11 — 21. Halos of this mass may host galax- 
ies of all types, not only the Sb galaxies consid- 
ered in this paper. One may argue that perhaps 
ellipticals, lenticulars, and Sa's form in high con- 
centration halos and low surface brightness (LSB) 
galaxies are hosted by low concentration halos. 
Thus, the range for "normal" Sb galaxies may be 
somewhat narrower than the full spread for all ha- 
los. Tentatively, we assume that the concentration 
should be in the range C = 10 — 17. 

The angular momentum of dark matter halos is 
also theoretically constrained. Halos in cosmologi- 
cal simulations have small amounts of angular mo- 
mentum, characterized by the dimensionless spin 
parameter 

A = J\E\V*/GM%? (18) 

where J is the angular momentum and E is the 
total energy of the halo. The spin parameter has 
a log-normal distribution, which does not depend 
(or has very little dependence) on cosmological pa- 
rameters, halo mass, or redshift (Barnes & Efs- 
tathiou 1987; Lemson & Kauffmann 1999; Vitvit- 
ska et al. 2001). We use the parameters given by 
Vitvitska et al. (2001). The maximum of the dis- 
tribution of A is at A = 0.035 and A has a 90% 
probability of being in the range 

A = 0.02 -0.10. (19) 

Just as in the case of the concentration, it is likely 
that there is a correlation of morphological type 
and spin, with Sb galaxies avoiding halos with 
both extremely low and high values of A. 

For each set of model assumptions, we can es- 
timate the spin parameter of its pre-collapse halo. 
After the fitting is done and we have the rotation 
curve, we find the angular momentum of the disk: 
■/disk =27r J V mt ^(r)r 2 dr, where £(r) is the sur- 
face density of the disk. We then assume that the 
specific angular momentum of the whole system 
j = J/M v i r is equal to the specific angular momen- 
tum of the present-day disk jdisk = Jdisk/A^disk- 
This gives us the total angular momentum J. Us- 
ing the parameters of the pre-collapse halo we es- 
timate the potential energy of the system and, as- 
suming virial equilibrium, find the spin parame- 
ter. The assumption of virial equilibrium should 



be quite reasonable because our systems are not 
interacting or merging. 

The assumption that j = jdisk is a sensible and 
conventional starting point, and should require lit- 
tle further discussion. We neglect the angular mo- 
mentum of the bulge in our calculations. This does 
not give large errors because the angular momen- 
tum of bulges in our models is much smaller that 
the angular momentum of disks. Because we re- 
strict our modelling to late type galaxies, bulge 
masses are also small (< 20% of the disk mass), so 
in the worst case we are making 20% error. The 
original angular momentum of the dark matter 
halo comes from gravitational (tidal) interactions 
with its environment. Thus, the dark matter and 
the gas experience the same torque in the process 
of halo assembly and should initially have (almost) 
the same specific angular momentum. The main 
uncertainty is what happens to the gas and its an- 
gular momentum from the moment it crosses the 
virial radius until it finally settles into the disk. As 
we shall find, a significant fraction of the gas must 
not end up in the disk. In our models, about half 
of the gas is not found in the disk or bulge. The 
models presented here do not address the question 
of where that gas is; whether it is still inside the 
virial radius, was expelled from the halo or was 
never in the halo. For simplicity, we assume that 
the specific angular momentum was not affected 
by whatever processes (presumably some form of 
feedback) were responsible for removing the gas 
from the disk or preventing its collapse. Another 
issue is the exchange of angular momentum be- 
tween the baryons and the dark matter. Some re- 
distribution of the angular momentum may have 
happened, but we do not know how significant it 
was. We consider two types of models, with and 
without exchange of angular momentum. In both 
cases, within a specific set of assumptions about 
the angular momentum exchange (or lack thereof), 
we can thus find the initial angular momentum of 
the gas which is presently in the disk. 

2.3. Baryonic compression 

We start with models with no exchange of an- 
gular momentum between different components. 
The effects of baryonic infall are treated follow- 
ing same approach as that of Blumenthal et al. 
(1986), Flores et al. (1993), and Mo, Mao, & 
White (1998). It is assumed that the compres- 
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sion occurs adiabatically and the angular momen- 
tum of each component is preserved as the baryons 
settle into the exponential disk. Assuming addi- 
tionally that the velocity anisotropy is preserved, 
one obtains the following relation: 



G [M b (r) + M dm (r)] r = GM^r^n 



(20) 



where j is the specific angular momentum, rj is 
the average radius of a dark matter particle before 
the baryonic compression, Mb(r) and M<j m are the 
baryonic mass and dark matter mass inside the 
final radius r and the Mh a \ {ri) is the halo mass 
before contraction, given by 

-Mhalol^iJ = 7Z • (21) 

* 'dm 

Equation (20) can be considered as an equation 
for the initial radius corresponding to a given 
final radius r. This equation is solved numerically. 

2.4. Exchange of angular momentum be- 
tween baryons and DM 

Models with exchange of angular momentum 
between the baryons and dark matter are more 
complicated. The exchange probably happens at 
late stages of the baryonic infall when the baryon 
density becomes large and a non-axisymmetric 
component may develop due to the excitation of 
spiral waves and/or bar- like modes. Giant molec- 
ular clouds may play some role in this process. 
Dynamical friction can then result in a transfer 
of angular momentum from the baryons to the 
dark matter. Because the dark matter gains angu- 
lar momentum, it moves further from the galactic 
center. Thus, the density of the dark matter in 
the central region decreases. In the early stages 
of galaxy formation, when most of the baryons 
were still in gaseous form, we might expect that 
those non-axisymmetric features were more preva- 
lent and more dynamically important then at the 
present time when most of the baryons are locked 
in stars. 

It is difficult to estimate the exact amount of 
angular momentum that would be lost by the 
baryons in such a situation, but we know that it 
cannot be very large. It is constrained by two 
factors: the angular momentum of the disk at 
present, and the condition that the initial spin pa- 
rameter should be close to the typical (median) 



value for halos. Models without the exchange of 
angular momentum give A « 0.02 — 0.03, which is 
already close to the most probable spin A = 0.035 
of a dark matter halo. If the disk loses some angu- 
lar momentum during its assembly, the estimate of 
the initial spin parameter increases almost linearly 
with the lost angular momentum. If too much an- 
gular momentum is lost, the spin parameter be- 
comes unacceptably large and we would be forced 
to conclude that our Galaxy formed from the ex- 
treme tail of the distribution. Still, because the 
distribution of spin parameter is quite broad, there 
is some room for angular momentum exchange. 
The spin parameter can be two times larger than 
for the no-exchange models and still be considered 
"typical" . Thus for our model with angular mo- 
mentum exchange, we assume that the disk loses a 
factor of 1.5-2 of its angular momentum during the 
collapse. For simplicity we assume that the forma- 
tion of the disk happens in two stages. During the 
first stage, when the baryons experience most of 
the collapse, they preserve their angular momen- 
tum. During the second stage the disk shrinks 
further, losing some of its angular momentum to 
the dark matter. 

We use the approach outlined in the previous 
section to compute the state of the system at the 
end of the first stage of adiabatic compression (in 
which angular momentum is conserved). We then 
consider a spherical shell of dark matter with ra- 
dius r, thickness dr 7 density pdm, and specific an- 
gular momentum 

j = rV c = y/G [M b (r) + M dm (r)]r. (22) 

It has a total mass dM dm , and total angular mo- 
mentum of 



dJ 
dM dm 



j dMdm 

4:wp dm r 2 dr . 



(23) 
(24) 



We estimate the amount of angular momentum 
lost by the baryons (disk and bulge) when a mass 
element of baryons AM& moves from radius r + dr 
to radius r. The baryons lose angular momentum 



dJ h = dM h 



(V c + ^dr)(r + dr) - V c r 



■ (25) 



This angular momentum is deposited into the dark 
matter, which produced the dynamical friction. 
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Thus, the final angular momentum of the DM shell 
is 

dJ f = dJ + dJ b . (26) 

Shells at different distances should be affected 
because the dynamical friction is not a local pro- 
cess. Nevertheless, for simplicity, we deposit all 
the angular momentum into the shell through 
which the disk or bulge clement moves. This is 
a reasonable approximation because (i) a signifi- 
cant fraction of the dynamical friction is due to 
the elements of dark matter closest to the sinking 
baryonic material, and (ii) the volume affected by 
the deposition of the angular momentum for all 
shells is very extended — up to 20 kpc. 

A shell of dark matter at radius r, which ac- 
quired specific angular momentum 



Jf 



dJ b /dM dm , 



(27) 



then moves to a new radius rf. Here jf is the final 
specific angular momentum and is related to the 
mass model by 



j) = GM(r f )r f . 



(28) 



Equating initial and final angular momenta for 
each shell, we get an implicit equation for the final 
radius rf. 



Jf = J 



1 + 



AAM 



A 
AM 



47r/9 dm r 3 
r^dVc 
V c dr ' 
M b , f - M b . 



1 + 



(29) 

(30) 
(31) 



Here M = Md m + Mb is the total mass inside 
a radius r. Eq. (29) is solved numerically. The 
solution also gives the mass inside a final radius r/. 
Eq. (29) has the same structure as eq. (20). The 
only difference is the term on the right-hand-side, 
which is the correction due to angular momentum 
deposition. 

We can get a rough estimate of the effect if 
we neglect the term A (which is close to unity) 
and introduce the average density excess produced 
by infalling baryons Ap = AM / (Anrf / 3) . If we 
further assume that M oc r (a good approximation 
for most radii), then eq. (29) takes the form: 



r f 



1 + 



Ap_ 

3pdm 



(32) 



ft is clear that during the initial stages of the col- 
lapse, when the density of the baryons was about 
ten times smaller than the density of the dark mat- 
ter, the exchange of angular momentum had little 
impact on the dark matter: rf « r^. The effect 
peaks at around Ap = 3pdm- At even larger val- 
ues of the density ratio, the approximation fails 
because a small amount of dark matter cannot ex- 
ert significant dynamical friction on large mass of 
baryons. Nevertheless, at the peak of its impor- 
tance, the effect is potentially quite large, with 
rf w 2ri, resulting in a decrease of the dark mat- 
ter density by a factor of ten. 

3. Observational Constraints 

3.1. Constraints for the Milky Way 

In this section we describe the observational 
data that we use to confront our models. 

Satellite dynamics and modeling of the Mag- 
ellanic Clouds (e.g. Zaritsky et al. 1989; Fich 
& Tremaine 1991; Lin, Jones, & Klemola 1995; 
Kochanek 1996) provide constraints on the mass 
of our galaxy on large scales (~ 50 — 100 kpc). 
The mass of the Milky Way inside 100 kpc is es- 
timated to be (5.5 ± 1) x 10 11 M from dynamics 
of the Magellanic Clouds, (Lin, Jones, & Klemola 
1995). Using constraints from the escape velocity 
and motions of satellite galaxies, Kochanek (1996) 
estimates the mass of the Galaxy inside 100 kpc 
to be (5 - 8) x lO n M . We adopt the following 
constraints on the mass inside lOOkpc (Dehnen & 
Binney 1998): 



M 



K<iookpc = ( 7±2 - 5 )x lollM o- (33) 



For the outer part of the rotation curve of our 
Galaxy we use the data summarized in Dehnen & 
Binney (1998). The radial circular velocity v r of 
an object at galactocentric coordinates I and b is 
related to the circular velocity at radius r, v c (r), 

by 



W(r/Ro) 



Re 



sin I cos b 



v c {r)-v c {R ), (34) 



where Rq is the Sun's distance to the galactic cen- 
ter. Data for H\\ regions (Brand & Blitz 1993) 
and for classical Cepheids (Pont et al. 1997) are 
used to derive W(r/Ro) for up to twice the solar 
distance. 
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The surface density of gas and stellar compo- 
nents at the solar radius is estimated by Kuijken 
& Gilmore (1989) to be: 

S st ars+ga S -48±8M Q p C - 2 (35) 

We also use another local constraint: the vertical 
force K z at 1.1 kpc above the galactic plane. Using 
data on the kinematics of K dwarfs, (Kuijken & 
Gilmore 1989, 1991) give the limit on the total 
density of matter inside 1.1 kpc: 

\K z (Ro, z = l.lkpc)|/27rG = (71 ± 6) Mgpc" 2 

(36) 

Oort's constants: 

A = (v c /R-dv c /dR)/2 (37) 
B = -(v c /R + dv c /dR)/2 (38) 

provide constraints on the circular velocity curve 
at the solar radius. We adopt the same values as 
Dchncn & Binncy (1998): 

A = 14.5 ±1.5 km sec" 1 kpc" 1 (39) 
B = -12.5 ±2 km sec" 1 kpc" 1 (40) 
A — B = 27 ± 1.5 km sec" 1 kpc" 1 . (41) 

The corresponding v c = (A — B)R n at the solar 
neighborhood is 200 — 240km/s. 

We use observations of terminal velocities 
^terminal to constrain the mass distribution and 
the rotational velocities inside the solar radius. If 
I is the galactic longitude, then for an axisymmet- 
ric model the terminal velocity is related to the 
circular velocity v c by 

^terminal = V C (Ro Sm I) - V c (Ro) Sm I, (42) 

where Rq is the distance to the galactic center. We 
make use of data from the HI surveys of Knapp, 
Stark, & Wilson (1985) and Kerr et al. (1986). 
The terminal velocities are compatible with the re- 
sults of Malhotra (1995). Dehnen & Binney (1998) 
present a more comprehensive comparison of dif- 
ferent surveys. 

To constrain the mass in the central parts of the 
MW, we use data on stellar motions. We make use 
of the results of Genzel et al. (2000), who found 
M = 3 x 10 7 M Q for the mass inside 10 pc and 
M = 10 7 M Q at 4 pc. Kinematics of OH/IR stars 
(Lindqvist, Habing, & Winnberg 1992) constrain 
the mass on the « 100 pc scale. 



More indirect constraints on the mass of 
baryons and dark matter in the central part of 
the galaxy can be obtained from the existence 
of persistent rapidly rotating bars and from the 
observed optical depth to microlensing events in 
the direction of the Galactic center. Based on 
the former, the work of Debattista & Sellwood 
(2000) suggests that dark matter must make up 
less than a quarter of the dynamical mass within 
one disk scale length. However, the halo models 
used in that work were rather unphysical n = 3 — 5 
polytropic cored models, and so it is difficult to 
know whether these results apply to cuspy CDM 
halos. From the microlensing constraints, Zhao & 
Mao (1996) concluded that all of the dynamical 
mass within 3.5 kpc must be in an optimal stel- 
lar bar to account for observed optical depth of 
(3 — 10) x 10" 6 . More recent data suggest that the 
optical depth may be much lower than this. Both 
of these pieces of evidence suggest that models 
with a minimal dark matter contribution in the 
inner part are to be preferred. We discuss how 
our models fare with respect to both of these con- 
straints in Section 5.2 and 6; however, we do not 
impose either of these constraints a priori. 

We use the K-band luminosity of our Galaxy 
derived by Drimmel & Spergel (2001): 

L K = 8.9 x 10 10 i© (43) 

corresponding to a magnitude of Mr- = —24.0. 
This value is consistent with the estimate of Mal- 
hotra et al. (1996), but it is larger than the value 
6 x 10 10 L Q adopted by Kent, D ame, & Fazio 
(1991). It is about 0.4-0.6 magnitudes fainter than 
the value implied by the K-band Tully-Fisher rela- 
tion of Tully & Pierce (2000), assuming a rotation 
velocity of 220 km/s for the Galaxy (see the dis- 
cussion in Section 4.3). 

3.2. M31: light and kinematics 

Throughout this paper, we assume that the dis- 
tance to M31 is 770 kpc, based on Cepheids and 
red clump stars (Freedman & Madore 1990; Ken- 
nicutt et al. 1998; Stanek & Garnavich 1998). In 
some earlier papers a smaller distance of 690 kpc 
was used. Results used below were rescaled to our 
adopted distance. 

The M31 galaxy has a more extended disk than 
our galaxy. Walterbos & Kennicutt (1988) studied 
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the surface brightness distribution of M31 in four 
colors, U, B, V, and R, tracing the disk to 20 kpc. 
The exponential length in the R-band was found 
to be 5.7 ±0.3 kpc, as compared with 2.5 — 3.5 kpc 
for our Galaxy. The scale length was found to be 
larger in bluer bands (6.4 kpc in B). We use the 
R-band value in our modeling because it is less af- 
fected by dust extinction and more likely to trace 
the overall stellar mass distribution. A significant 
fraction of the light (1/3 to 1/2, depending on 
the band) is from the bulge. The bulge of M31 
is triaxial, as evidenced by a significant twist in 
isophotes and complex gas velocity patterns. Like 
the Milky Way, M31 has a nucleus as well. In fact, 
the nucleus can be resolved into two central peaks 
separated by about 1 pc. A black hole of mass 
<~ 3.7 x 1O 7 M is harbored in the less luminous 
peak (P2). 

We model the bulge of M31 with a scaled- up 
version of the Milky Way bulge (eqs. 1-4) and disk, 
using a disk scale radius of = 5.7kpc. Wc in- 
clude a black hole with mass 3.7 x 10 7 M — more 
than ten times more massive than the Milky Way 
black hole. R-band photometry along the major 
axis of M31 is taken from Walterbos & Kennicutt 
(1987). We also use r— band photometry from 
Kent (1987). The latter results are shifted up 
by 0.45 mag to match the R-band surface bright- 
ness. We apply a correction of 0.25 mag to both 
the bulge and disk to account for extinction in 
our Galaxy, and apply an additional correction of 
0.74 mag to the disk for internal extinction (Kent 
1989). 

We use the CO observations of Loinard, Allen, 
& Lequeux (1995) to determine the rotational ve- 
locity inside the central 10 kpc. HI measurements 
extend the rotation curve to 30 kpc (Brinks & Bur- 
ton 1984). The CO and HI velocities agree in the 
overlapping central region of the galaxy. The cen- 
tral 2 kpc region of M31 shows anomalously high 
velocities and a misalignment of the major axes 
of the bulge and disk. Those anomalies are at- 
tributed to the presence of a triaxial bar/bulge 
in M31 (e.g. Stark 1977; Stark & Binney 1994; 
Berman 2001). Just as in the case of the Milky 
Way, we use a simplified spherical model of the 
central bulge. The model can not reproduce the 
details of the circular velocity function in the cen- 
tral 2-5 kpc, but it gives useful constraints on the 
mass models. 



Additional constraints are obtained from stel- 
lar kinematics in the central region of M31. We 
use the estimates of the rotation velocity v mt and 
the line-of-sight velocity dispersion a from Kor- 
mendy & Bender (1999, Table 2). We then esti- 
mate the mass, assuming an isothermal distribu- 
tion with isotropic velocities M = 2((j 2 + v^ ot )r/G. 
At distances 25"-50" (w 100 - 200 pc) from the 
center of M31 the rotational velocity is small, the 
velocity dispersion is almost constant (150km/s- 
160km/s) and the assumed density profile of the 
nucleus is close to r~ 2 . Thus, our approximation 
seems to be quite reasonable for these radii. 

4. Models of the Milky Way and M31 
4.1. Comparison with Observations 

Our goal is to find whether a set of models, 
which we consider plausible or theoretically well- 
motivated, are compatible with the observational 
constraints. We do not attempt to find best-fit 
models, but present a range of models and dis- 
cuss whether they satisfy the constraints we have 
adopted. Once we find by trial and error that a 
model is compatible with the constraints, we do 
not attempt to further improve the quality of the 
fit by fine-tuning the parameters of the model. 

Tables 1-3 list the parameters of our models. 
Most of the parameters and values have already 
been introduced and discussed in sections 2 and 
3. The last row in each table gives the ratio of 
the sum of disk and bulge masses to the mass of 
baryons expected inside the virial radius (/fcM v j r ). 
We assume the universal ratio of baryons to dark 
matter is fb = flb/(fldm + ^fc) = 0.1. Our results 
can be easily adjusted to any other ratio. 

With the exception of models A 2 and B 2l all 
the models pass the constraints we have adopted. 
Typically the fitting was done by first assuming a 
mass and concentration for the halo and then by 
tuning the disk parameters. 

In Models A 2 and B 2 we used a different ap- 
proach. Here we started by assuming a disk mass 
of M 3 = 6 x 10 10 M Q . This mass is only 20-50% 
more massive than the value used in our favored 
models Ai and B\. It is often quoted as the fidu- 
cial mass of the disk (Binney & Tremaine 1987). 
We then tried to find an acceptable fit to the data 
by changing parameters of the halo. Even assum- 
ing an extremely low halo mass and unrealistically 
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Table 1 

The Milky Way galaxy: No exchange of angular momentum 



Parameter 


Constraints 


Model A x 
Favored 


Model A 2 
Max. Disk 


Model A 3 
Max. C 


Model A± 
Max. halo 


Virial mass M V - 1T (M ) 


- 


1.0 x 10 12 


0.71 x 10 12 


1.0 x 10 12 


2.0 x 10 12 


Virial radius r v i r (kpc) 


— 


258 


230 


258 


325 


Halo concentration C 


10-17 


12 


5 


17 


10 


Disk mass M 3 (cq. (3); M Q ) 




4 x 10 10 


6 x 10 10 


3.5 x 10 10 


4.0 x 10 10 






llSyl n 10 


i 9 v i n 10 


n 7 v i n 10 


nsyi n 10 

U.O A _LU 


Stellar mass r < 3.5 kpc (M©) 




1.9 x 10 10 


3.1 x 10 10 


1.5 x 10 10 


1.9 x 10 10 


Disk exponential scale length 


2.5-3.5 


3.5 


3.0 


3.5 


3.5 


rd (kpc) 












Maximum circular velocity 




228 


246 


235 


237 


Vc.max (km/s) 












Max. halo circular velocity 




163 


123 


178 


197 


no compression (km/s) 












Solar distance, Rq (kpc) 


7-8.5 


8.0 


8.5 


8.5 


8.5 


Baryon surface density at Rq, 


48 ±8 


53 


62 


40 


46 


(M pc- 2 ) 












Total surf, density within 1.1 kpc 


71 ±6 


75 


74 


63 


69 


at R G (M Q pc~ 2 ) 












Oort's constants A — B 


27 ± 1.5 


26.8 


28.1 


27.6 


27.9 


(km s _1 kpc^ 1 ) 












Mass inside 100 kpc (10 11 M Q ) 


7.5 ±2.5 


5.8 


3.8 


6.1 


9.0 


Disk+bulge mass-to-light ratio 




0.54 


0.81 


0.47 


0.54 


(M/L) K (M /L ) 












Spin parameter A 


0.02-0.10 


0.031 


0.022 


0.037 


0.018 


M dm /M disk+ buigo(r < 3kpc) 




0.95 


0.40 


1.31 


1.02 


Fraction of "galactic" baryons 




0.48 


1.00 


0.42 


0.24 
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Table 2 

The Milky Way galaxy: Models with the exchange of angular momentum 



Parameter 


Constraints 


Model ri\ 


Model D2 


Model ±>3 


Model ±>4 






Favored 


iviax. uisk 


Large Mass 


oman rLxcn. 


Virial mass iVf v ir (Mq) 




1.0 x 10 12 


0.71 x 10 12 


1.5 x 10 12 


1.0 x 10 12 


Virial radius r V j r (kpc) 




258 


230 


295 


258 


Halo concentration C 


10-17 


12 


10 


10 


12 


Disk mass M 3 (M ra ) 




5 x 10 10 


6 x 10 10 


5 x 10 10 


5 x 10 10 


Mass Mi 4- Mo (MrO 




1 x 10 10 


1.2 x 10 10 


1 x 10 10 


1 x 10 10 


Stellar mass r < 3.5 kpc (M ) 




2.7 x 10 10 


2.8 x 10 10 


2.7 x 10 10 


2.6 x 10 10 


Disk exponential scale length 


2.5-3.5 


3.0 


3.5 


3.0 


3.0 


Vri (kpc) 












Intermediate exponential scale 


- 


6.0 


7.0 


6.0 


4.5 


rd (kpc) 












Maximum circular velocity 




223 


216 


225 


234 


Vc.max (km/s) 












Max. halo circular velocity 




163 


139 


178 


163 


no infall (km/s) 












Solar distance Rq (kpc), 


7-8.5 


8.5 


8.5 


8.5 


8.5 


Baryon surface density at Rq, 


48 ±8 


52 


69 


52 


52 


(M pc- 2 ) 












Total surf, density within 1.1 kpc 


71 ±6 


68 


79 


70 


72 


at R & (M pc~ 2 ) 












Oort's constants A — B, 


27 ± 1.5 


26.5 


25.4 


26.5 


28.0 


(km s _1 kpc -1 ) 












Mass inside 100 kpc , 


7.5 ±2.5 


6.0 


4.6 


7.7 


6.0 


(lO n M ) 












Bulge+disk mass-to-light ratio 




0.67 


0.81 


0.67 


0.67 


(M/L) K (M /L ) 












Spin parameter A 


0.02-0.10 


0.058 


0.076 


0.044 


0.050 


M dm /M disk+bu i ge (r < 3kpc) 




0.24 


0.14 


0.24 


0.31 


Fraction of "galactic" baryons 




0.6 


1.0 


0.4 


0.6 
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Table 3 

Parameters of Models for the M31 galaxy 



Parameter 


Model Ci 


Model C 2 




No Exchange 


With Exchange 



Virial mass M v ; r (M©) 


1.60 x 10 12 


1.43 x 10 1 


Virial radius r v ; r (kpc) 


300 


290 


Halo concentration C 


12 


12 


Disk mass M d (M Q ) 


7.0 x 10 10 


9.0 x 10 10 


Bulge mass M b (M Q ) 


1.9 x 10 10 


2.4 x 10 10 


Disk exponential scale r d (kpc) 


5.7 


5.7 (8.55) 


Stellar mass r < 3.5 kpc (M ) 


2.8 x 10 10 


3.6 x 10 10 


DM mass r < 3.5 kpc (M Q ) 


2.5 x 10 10 


1.8 x 10 10 


Bulge (M/L) R (M /L ) 


3.0 


3.8 


Disk (M/L) R (Mq/Lq) 


0.93 


1.2 


Maximum circular velocity VJ; )max (km/s) 


269 


262 


Max. halo circular velocity Vhaio,max 


184 


176 


without compression (km/s) 






Mass inside 100 kpc radius (M Q ) 


8.5 x 10 11 


8.1 x 10 11 


Spin parameter A 


0.036 


0.057 


Fraction of "galactic" baryons 


0.56 


0.80 



low concentration, both models fail. In Model 
A2, the maximum circular velocity and the sur- 
face density at the solar radius are both too large. 
Model B 2 is a bit better, but still the local surface 
density is too large, the angular momentum is un- 
comfortably large, and the mass inside 100 kpc is 
too small. We can decrease the local surface den- 
sity by reducing the exponential length of the disk. 
But that has the side effect of increasing the disk 
mass inside the central 5 kpc region. As a result, 
the rotation velocity gets unacceptably large and 
the model fails again. 

Some cautions should be given regarding the 
bulge mass and the mass distribution inside the 
central «3 kpc. The values Mi, M 2 and M 3 are 
the masses of three inter-penetrating components, 
which are the dominant components in the regions 
of the nucleus, the bulge/bar and the disk respec- 
tively. The p 3 component is a constant-height thin 
exponential disk, which remains as such even in- 
side the central «3 kpc region occupied by the 
galactic bar. This is unphysical and is just a sim- 
plification for the convenience of the fitting. In 



reality the disk should be mixed with the bar- 
bulge in the central region and should constitute a 
single component with the distribution of the cen- 
tral non-axisymmetric bar. These complexities are 
only important for a transitional region between 
the bulge and the disk at radii 2-4 kpc, where the 
disk mass is comparable to the bulge mass. Inside 
2 kpc the disk is much smaller than the bulge and 
it does not matter how we treat the disk. 

Some important quantities for the best mod- 
els of the Milky Way and M31 are presented in 
Figures 1-5. Figure 1 compares the observed 
terminal velocities in our Galaxy with the predic- 
tions of model A\ . The quality of the fit decreases 
at small angles because of non-circular motions 
produced by the central bar. Even larger dis- 
tortions are found in the rotation curve of M31, 
which has a large bar. The usual way to deal 
with the deviations is to ignore the data at small 
angles. For example, Dchncn & Binncy (1998) 
ignore data for angles smaller than 17°, which 
corresponds to « 2.5 kpc. Oiling & Merrifield 
(1998) used the data of Malhotra (1995), which 
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Fig. 1 . — Dependence of terminal velocities on galac- 
tic longitude I. The full curve is for model A\. Sym- 
bols show observational data from the HI measure- 
ments of Knapp et al. (1985) (circles) and Kerr et 
al.(1986) (triangles). At small angles the deviations 
from circular velocities are expected to be large due to 
the central bar. This is clearly seen at I < 20°. 

also start at w 2.5 kpc. Without a more realistic 
(non-axisymmetric) model of the bar one cannot 
do better than this. The problem is that we can 
not exclude the possibility that small 5 — 10 km/s 
disturbances related to the bar can be felt even 
at larger radii. This could bias our parameter es- 
timates, which are sensitive to what happens at 
those radii. For this reason, we still consider the fit 
produced by model A\ to be acceptable in spite of 
the fact that for longitudes 20°— 30° the model val- 
ues exceed the observational points by 5 — 10 km/s. 

Figure 2 shows the rotation curve of our Galaxy. 
The observational data are compared with mod- 
els A\ and B\ . We also show the contributions of 
different components. For both models the dark 
matter dominates in the outer parts of the Galaxy, 
but the radius at which the contribution of the 
dark matter equals that of the baryonic compo- 
nent are very different. For model A\, the masses 
contributed by baryons and dark matter are equal 
at 3.5 kpc. For model B\ this radius is 11 kpc. The 
difference is due to a combination of two factors: 
expulsion of dark matter by dynamical friction and 



an increase in the assumed disk+bulge mass. In 
the central 5 kpc of the Galaxy the contribution 
of the dark matter relative to the disk and bulge 
is very different for models A\ and B\ . In the case 
of model A\ and in all other models with no ex- 
change of angular momentum, the dark matter is 
never a strongly dominant component, but it al- 
ways contributes 40%-60% of the total mass. In 
other words, these models have sub-maximal disks 
(Bottema 1997; Courteau & Rix 1999). In the 
case of models with exchange of the angular mo- 
mentum, there is little dark matter in the central 
region of our Galaxy. 

Figure 3 shows the distribution of mass in the 
Milky Way galaxy for a very large range of scales, 
from 5 pc to 200 kpc. The following observa- 
tional constraints are used. The first two points 
at 5-10 pc are from studies of stellar radial ve- 
locities and proper motions in the galactic cen- 
ter. Squares are based on kinematics of OH/IR 
stars (Lindqvist, Habing, & Winnberg 1992). The 
point at 3.5 kpc is based on the Zhao (1996) model 
of the bar. Because the model was compared 
with the data on stellar kinematics (inner rota- 
tion curve and radial velocity dispersion), it gives 
a constraint on the total mass: 4 x 1O 1O M , with 
an uncertainty of about 20%. For the next data 
point at 8.5 kpc we simply assume that the cir- 
cular velocity is 220 ± 20 km/s, which covers the 
whole range of reasonable values. We then esti- 
mate the mass as M = v 2 r/G. The last observa- 
tional point is the constraint from the motions of 
satellite galaxies discussed in section 3. The cen- 
tral data points were not used in either our fitting 
or in the analysis of the bulge (Zhao 1996). Never- 
theless, they come fairly close to the extrapolation 
of our model into the very center of our Galaxy. 
The theoretical curves for our favored models A\ 
and B\ are very close to each other, which is not 
surprising because they fit the same data and have 
the same global dark matter content. The largest 
deviation of the models from the data is for the 
mass inside 100 pc, where the observational es- 
timate is twice larger than the prediction of the 
models. Even at this point the disagreement is 
not alarming because the observational data are 
likely more uncertain than the formal error. 

What is remarkable about Figure 3 is that it 
spans more than 5 orders of magnitude in radius 
and mass. It is encouraging that, without fine- 
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R(kpc) 

Fig. 2. — Rotation curve for our favorite models Ai 
(no exchange of angular momentum) and B\ (with the 
exchange) . Note that the dark matter dominates only 
in the outer part of the Milky Way. Symbols show 
observational data from HI measurements of Knapp et 
al. (1985) (circles) and Kerr et al.(1986) (triangles). 



tuning, our models are consistent with observa- 
tions of the dynamical mass of the MW over this 
huge range. 

Finding an acceptable model for M31 was rel- 
atively easy because there are much less data. In 
particular, we do not have kinematic constraints 
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Fig. 3. — Mass distribution of the MW galaxy for 
Model Ai (full curve) and model Bi (dashed curve). 
The large dots with error bars are observational con- 
straints. From small to large radii the constraints are 
based on: stellar radial velocities and proper motions 
in the galactic center; radial velocities of OH/IR stars; 
modeling of the bar using DIRBE and stellar veloci- 
ties; rotational velocity at the solar radius; dynamics 
of satellites. 



for the disk, which would be equivalent to con- 
straints at the solar position in our Galaxy. Our 
model seems to reproduce reasonably well the dy- 
namical mass of M31 from 100 pc to «100 kpc. 
Our model does not produce the very large wig- 
gles exhibited by the observed rotation curve. The 
wiggles at 5 kpc and 9 kpc are likely due to non- 
circular motions induced by the bar and, thus, as 
discussed before, cannot be reproduced by any ax- 
isymmetric model. The bulge of M31 is almost 
twice as massive as the bulge of our Galaxy. It 
is also slightly (30%) more compact. The disk 
of M31 is also more massive, but it is more ex- 
tended. As a result, in the central 5 kpc of the 
M31 the bulge is a much more dominant compo- 
nent as compared with the bulge of our Galaxy. 

The surface brightness profile in the R-band, 
shown in figure 5, is used as an additional con- 
straint. An accurate fit (the same as for the mass 
modeling) is obtained for stellar mass-to-light ra- 
tios of M/L = O.93M O /L and M/L = 3M /L 
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for the disk and the bulge respectively. These re- 
sults are quite consistent with the expectations 
from stellar population synthesis models for a 
galaxy with the B — R color of M31, which imply 
an overall stellar mass-to-light ratio (M/L)r = 
O.85M Q /L (Bell & de Jong 2001). The bulge 
M/L is likely overestimated, and would be re- 
duced by a correction for internal dust absorption. 
For example, internal absorption of 0.25 mag in- 
side the bulge would reduce M/L to 2.4. At a first 
glance our estimates of the M/L ratios are signif- 
icantly different from those of Kent (1989), who 
give M/L = 5M /L and M/L = 1O.4M /L Q 
for the bulge and the disk. Most of the differ- 
ences come from differences in bands (r instead 
of R) and from the fact that Kent's luminosities 
were not corrected for absorption. If we include 
those corrections (see Section 3.2) and also scale 
the results to the distance of 770 pc, Kent's es- 
timates become M/L = 2.5M Q /L0 for the disk 
(M/L « 6.7 in B) and M/L = 2.4M Q /L for the 
bulge. Such a large M/L ratio for the disk, more 
typical of an old bulge population, seems a bit 
problematic. The main reason why Kent (1989) 
has such a large M/L lies in the assumed profile 
of the dark matter: constant density through the 
whole galaxy. This results in the extreme case of 
a maximal disk, and very little dark matter inside 
20-25 kpc. 

4.2. Effects of compression by baryonic in- 
fall 

The sinking of baryons into the central part of 
a galaxy increases the depth of the gravitational 
potential. This leads to an increase of the dark 
matter density. To some degree, this process pro- 
vides coupling of the baryons and the dark matter. 
As a result, in the case of adiabatic compression 
without the exchange of the angular momentum, 
it appears difficult to find a case where one of the 
components (dark matter or baryons) is signifi- 
cantly larger than the other. Figure 2, top panel, 
gives an example of this effect. From 100 pc to 
about 6 kpc the contributions of baryons and dark 
matter to the circular velocity are about equal. 

We use two simple models to study the effect of 
adiabatic compression on the ratio of dark matter 
and baryons. In the first case we take Model Ai 
and change the mass of baryons while keeping the 
disk scale length and the bulge-to-disk ratio con- 
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Fig. 4. — Mass distribution (top panel) and rotation 
velocity (bottom panel) of the M31 galaxy for Model 
Ci. The large dots with error bars are observational 
constraints. In the bottom panel the circles show re- 
sults of CO (r < 10 kpc) and HI (r > 10 kpc) obser- 
vations. The large dip at around 5 kpc is likely due 
to non-circular motions induced by the bar. Obser- 
vational data points in the top panel at small radii 
are from stellar motions in the nucleus. Vertical bars 
correspond to 20% errors in mass. Data points at 
7, 15, and 30 kpc correspond to circular velocities of 
240 ± 10, 270 ± 20, and 240 ± lOkm/s. 



stant. In the second case, we vary the scale length 
of the disk, keeping the mass of the disk the same 
as in Model Ai; we ignore the bulge to make the 
case more clear. Figure 6 shows the ratio of the 
circular velocity produced by the dark matter to 
that of the baryons. 

Figure 6 illustrates that for a realistic range of 
parameters the contribution of baryons can not be 
significantly larger that of the dark matter. Even 
in the case when disk is four times more massive 
or four times more compact than in Model A±, 
the contribution of the dark matter to the rota- 
tion curve was at least 1/3. In other words, in hi- 
erarchical models with adiabatic compression and 
no exchange between components it is difficult to 
have a maximum disk. 

It is interesting to note that the opposite is also 
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Fig. 5. — Surface brightness of M31 in the R-band 
on linear (bottom panel) and logarithmic (top panel) 
scales. Deviations from the observational results are 
less than 0.2 mag. 

true: in the central part of a galaxy it is difficult 
to have a model where the mass of the dark matter 
is much larger than that of baryons. The model 
with a four times more extended disk in the top 
panel (x4) is an example of what should represent 
a low surface brightness galaxy. Its surface density 
is 16 times smaller than for our Galaxy, and its 
mass in stars is only 1% of the virial mass, yet the 
baryons at the solar radius contribute about 1/4 
of the rotational velocity 

In view of this self-adjusting behavior of the 
dark matter, it is easier to understand the profiles 
of different components in realistic models. The 
top panel in Figure 7 shows the dark matter and 
the combined disk and bulge density profiles. For 
comparison we also show the dark matter profile 
before the baryonic infall. Note that the change 
in the dark matter profile is quite large, yet the fi- 
nal profile manages to follow the baryonic density 
very closely in the central region. The bottom 
panel shows different components for model Bi. 
Each component has a complicated shape. Dif- 
ferent components become prominent at different 
radii. The net effect is rather remarkable — the 
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Fig. 6. — The ratio of the circular velocity produced 
by the dark matter to that of the baryons for two mod- 
els of adiabatic compression. For the top panel the ex- 
ponential disk scale length rd was changed by a factor 
indicated in the plot. The factor of unity (xl) is for 
the assumption of a disk mass 4 x 10 10 M Q and scale 
length Td = 3.5 kpc. This model has no bulge. Curves 
in the bottom panel show the same quantity when the 
disk scale length is held fixed at rd = 3.5 kpc the total 
baryonic mass is varied. Parameters of the xl curve 
are the same as for Model A\ . 

sum of all components is very close to a power-law 
with slope -2. The deviations from this power law 
are less than 10% on scales of 1-20 kpc. One may 
argue that this occurs because we are fitting a ro- 
tation curve that is nearly flat, and so we get a r~ 2 
density profile. This is not exactly true. The ob- 
served deviations are quite substantial. The final 
fit is much closer to a pure power law than what 
is required by the observations. This "disk- halo 
conspiracy" is the natural result of the coupling of 
baryons and dark matter via the adiabatic com- 
pression. 
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Fig. 7. — Density distribution for different compo- 
nents in models A\ (top panel ) and Bi (bottom panel. 
Different curves show: the total density (full) and dark 
matter (short dash). Top: The total bulge+disk den- 
sity is shown by the dot-dashed curve. Inside w 3 kpc 
the baryons and the dark matter closely follow one 
another. This illustrates the self-adjusting nature of 
the compression. The dotted curve shows the dark 
matter before the compression by baryons. The com- 
pression is clearly a dominant effect inside the central 
10 kpc. Bottom: The bulge is shown by the long- 
dashed curve; the disk is shown by dot-dashed curve. 
The components have very different profiles, but their 
sum is very close to a power-law with slope -2, shown 
by the straight line. This "disk-halo conspiracy" is the 
natural result of the adiabatic dark matter compres- 
sion. 



4.3. Constraints from the Global Lumi- 
nosity Function 

If the luminosity of our Galaxy is typical of 
that of other disk-type galaxies that form in halos 
of the same mass, this also has implications for 
the global luminosity function. We now investi- 
gate the number density of galaxies like the Milky 
Way, using the revised Press-Schechter approxi- 
mation of Sheth & Tormen (1999). The Sheth- 
Tormcn approximation provides us with the num- 
ber density of dark matter halos as a function of 
their mass. Once we have pinned down the lu- 
minosity of the Milky Way, we can calculate the 
mass-to-light ratio within the virial radius of the 
dark matter halo hosting the galaxy in a given 
model. We can then convert the halo number den- 
sity per mass interval to a galaxy number density 
per magnitude interval, which may be compared 
with the observed luminosity function. We neglect 
the contribution from galaxies that are "satellites" 
in much larger mass halos (i.e., cluster galaxies), 
and assume that the adopted mass-to-light ratio 
holds for disk galaxies only. We then compare 
this prediction with the K-band luminosity func- 
tion for all galaxies and for late-type galaxies from 
the 2MASS survey (Kochanek et al. 2001). 

We use two different approaches to estimate the 
K-band luminosity of the Galaxy. The total lu- 
minosity of our Galaxy is notoriously difficult to 
measure directly because we live in the middle of 
it, and obscuration is important even in the K- 
band. Drimmel & Spergel (2001) recently per- 
formed a detailed analysis in which they estimated 
the K-band luminosity of the Milky Way by fitting 
a model to the COBE DIRBE photometric data. 
They found Mk = —24.0 for the total luminosity 
of the Galaxy. 

The other approach is to assume that the 
Galaxy is "typical" and should lie on the Tully- 
Fisher relation (TFR). Unfortunately, this ap- 
proach is also uncertain because there is little 
Tully-Fisher data in the K-band and the results in 
the literature are somewhat inconsistent. If we use 
the TFR of Malhotra et al. (1996), assuming that 
the velocity width of the Milky Way AV = 2V C = 
440 km/s, then we obtain Mk = —23.37, about 
0.7 magnitudes fainter than the direct DIRBE 
measurement. If we use instead AV = 480 
km/s, which Malhotra et al. obtain by scaling 
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from M31 (AV(MW) = (220/250) AV(M31), with 
log 10 (AV(M31)) = 2.737), then we get -23.69, 
only 0.3 magnitudes fainter than the direct mea- 
surement. This would place the MW about l-cr 
brightwards of the TFR, which seems not unrea- 
sonable. However, the K'-band TFR recently pub- 
lished by Tully & Pierce (2000) and Rothberg et 
al. (2000) gives M K = -24.4 for AV = 440 km/s. 
The slope is very similar to the relation of Malho- 
tra et al., but the zero-point is almost one full mag- 
nitude brighter. This problem is only exacerbated 
if we consider that the Cepheid calibration data 
used by Malhotra et al. is consistent with a Hub- 
ble parameter of H a = 71 km/s/Mpc, while the 
ground-based K' data imply H = 81 km/s/Mpc. 
If we scale the zero-point of the ground-based data 
to the lower value of Ho ~ 70 km/s/Mpc, all the 
galaxies are further away and therefore are even 
brighter, making the discrepancy worse (the Milky 
Way would then have M K = -24.7). 

Malhotra et al. (1996) base their TFR on pho- 
tometry from DIRBE and have a total of seven 
galaxies in their sample (including the Milky 
Way). The ground-based K'-band sample consid- 
ered by Rothberg et al. (2000) and Tully & Pierce 
(2000) are from the Ursa Major sample (see also 
Verheijen 2001) and additional galaxies from the 
Pisces filament, and consists of a total of 69 galax- 
ies with 4 galaxies having Cepheid distances. All 
of these galaxies have multi-band photometry in 
B, R, I and K' and the TF relations in the optical 
bands agree well with those from other samples, 
including the very well-studied I-band. Therefore 
it seems that the zero-point for the ground-based 
sample should be more reliable. This suggests 
that perhaps there is some sort of calibration off- 
set in the DIRBE photometry, and perhaps calls 
into question the validity of using even the Drim- 
mel & Spergel (2001) result. We shall consider 
the implications of using three different choices 
for the total K-band luminosity of the Galaxy: 
the direct DIRBE result, the raw ground-based 
TFR result, and the ground-based TFR result 
scaled to H = 70 km/s/Mpc. The three sets of 
squares plotted in Figure 8 correspond to these 
three values. 

In each case the higher of the two connected 
squares shows the number density that we would 
obtain if every halo hosted a disk galaxy. The 
lower point shows the value obtained if only half of 



the halos host disk-type galaxies (probably a rea- 
sonable lower limit, based on semi-analytic mod- 
elling). The points can easily be scaled up or 
down for different assumptions about the fraction 
of disk-hosting halos, but for reasonable assump- 
tions the result should lie between the lower and 
upper points. Unfortunately, the interpretation of 
the luminosity function constraint is quite sensi- 
tive to the uncertain value of the MWs total lu- 
minosity. For the bright normalization (from the 
TFR of Rothberg et al. (2000) scaled to H = 70 
km/s/Mpc), the number density of halos harbor- 
ing MW galaxies is about a factor of two higher 
than the observations. For the faint normalization 
(the direct value from Drimmel & Spergel 2001), 
the number density of MW halos is marginally 
consistent with the observed luminosity function. 

It is of some concern that for the Milky Way 
normalization that we regard as the most reliable 
(the ground-based TFR scaled to the standard 
value of Ho), the luminosity function constraint 
is violated. If this normalization turns out to be 
correct, then this may suggest that the Milky Way 
is harbored by a halo with a larger virial mass. 
The number density of dark matter halos drops 
fairly rapidly with virial mass so this relaxes the 
constraint somewhat. Figure 9 shows the number 
density of dark matter halos as a function of their 
virial mass, normalized by the observed number 
density of galaxies at the assumed luminosity of 
the Milky Way from Kochanek et al. (2001). We 
see that the faint and even the intermediate Milky 
Way luminosity can be accomodated fairly com- 
fortably within the mass range allowed by our dy- 
namical modeling (10 12 M < M vir < 2xl0 12 M Q ; 
with brighter Milky Ways favoring larger virial 
masses), however, the bright Milky Way normal- 
ization is too high even at the upper virial mass 
limit of 2 x 10 12 M Q . 

5. Discussion 

There are numerous aspects of these models of 
the MW and the M31 galaxies which are impor- 
tant to consider. We start with very small scales 
and proceed to larger ones. 

5.1. Central region and the black hole 

It is now quite well established that black holes 
of mass ~ 10 6 — 10 8 M Q reside in the centers 
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Fig. 8. — The K-band galaxy luminosity function. 
The curved lines show the observed K-band luminos- 
ity function from the 2MASS survey (Kochanek et 
al. 2001) scaled to H = 70 km/s/Mpc; the solid 
line shows the results for all types of galaxies and 
the dashed line shows the results for late-type galaxies 
only. Points show the prediction of the number den- 
sity of halos hosting "Milky Way" galaxies, using the 
halo virial mass from model Ai or Bi and an assumed 
Milky Way luminosity. Solid symbols show the value 
obtained if all halos are assumed to host disk galax- 
ies, and are to be compared with the solid lines. Open 
symbols show the result obtained if only half of the ha- 
los host disk galaxies, and are to be compared with the 
dashed lines. The three sets of connected points are 
for three different assumed values for the total K-band 
luminosity of the Milky Way (see text). 

of most galaxies. The case for a 2.6 x 10 6 M Q 
black hole at the center of the Milky Way is de- 
manded by rising radial velocity and proper mo- 
tion dispersions within a few pc of the center 
(Ghez et al. (2000) and references therein). Re- 
cently, Gondolo & Silk (1999) showed that if the 
BH is grown adiabatically, the central dark matter 
will be drawn into a dense spike, with density ex- 
ceeding 10 s M Q /pc 3 . This would make the Galac- 
tic center very luminous in terms of neutrino flux 
if the dark matter is made of annihilating neutrali- 
nos. Ullio, Zhao, & Kamionkowski (2001) however 
argue that it is unlikely for the BH to grow all its 



Fig. 9. — Constraints on the virial mass of the Milky 
Way from the observed luminosity function. The top 
panel shows the ratio of the predicted to observed 
number densities of Milky Way galaxies as a function 
of the assumed average virial mass of dark matter ha- 
los hosting galaxies with the luminosity of the Milky 
Way. The three lines correspond to the three values of 
the Milky Way luminosity discussed in the text; where 
the dotted line is the bright normalization, the solid 
line is the intermediate normalization and the dashed 
line is the faint normalization. The plotted line should 
be less than unity in order for the model to satisfy the 
luminosity function constraint. Larger masses have 
lower number densities and therefore more easily sat- 
isfy the constraint. The bottom panel shows the halo 
mass-to-light ratio (i.e. the halo virial mass divided by 
the total galaxy luminosity) as a function of halo mass. 
Again, the three lines correspond to the three different 
Milky Way luminosities considered in the text. 

mass while at rest in the center because of galaxy 
merging in CDM models. A spiraling-in BH would 
actually "heat" and reduce the central density of 
the dark matter. The seed is likely to be mas- 
sive (> 1O 5 M ), because otherwise it would not 
have enough time to spiral into the center by dy- 
namical friction. Realistic models of loss-cone cap- 
ture show that the rate of stellar capture is around 
10 6 M per Hubble time (Ma gorrian & Tremaine 
1999). This implies that the BH could at most 
double its mass during the adiabatic phase. 
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Following Ullio, Zhao, & Kamionkowski (2001) 
we take into account the effect of adiabatic growth 
from a seed BH of half the final value. The flux is 
very likely to be small because in our models the 
final peak of the dark matter density outside the 
event horizon of the BH is well below 10 8 M Q /pc 3 , 
which is the amount needed for significant annihi- 
lation. 

5.2. Microlensing Constraints 

Microlensing events toward the Galactic cen- 
ter provide a lower limit on the mass of baryons 
within the inner part of the Galaxy. Combined 
with the observed terminal velocities, this places 
upper limits on the mass of dark matter that can 
be present and implies that the inner Galaxy must 
be mostly baryonic with a massive bulge and disk. 
The microlensing events toward the Galactic bulge 
are still not completely understood observationally 
and theoretically, in particular the nature of two 
long-duration events near the I = 3°, b = —3° 
field. Several values for the spatially averaged op- 
tical depth are reported with large discrepancies 
in the mean values. The earlier values are too 
high for axisymmetric bulge and disk models, and 
are barely consistent even with a massive bar. A 
2.2 x 10 10 M bar pointing about 13° - 20° away 
from the Galactic center-Sun line plus a full disk 
with a mass of 1.7 x 10 10 M Q inside 3.5 kpc are re- 
quired by detailed dynamical models of the COBE 
map and by the kinematics of the bulge stars 
(Zhao 1996). Zhao, Spergel, & Rich (1995) found 
that the bulk of the lenses are on the near side of 
this elongated bar. This bar plus disk model can 
produce an optical depth consistent with the lower 
la error bar of the observed microlensing optical 
depth of (3.3 ± 1.2) x 10~ 6 for all stars observed by 
OGLE (Udalski et al. 1994) toward Baade's win- 
dow (I = 1° and b = -4°). Zhao, Rich, & Spergel 
(1996) found that this model with a flat IMF with 
very few brown dwarfs can fit the microlensing 
event duration distribution. Peale (1998) stud- 
ied the optical depth and the distribution of dura- 
tions of lensing events for the bulge with the shape 
used in our paper. Both statistics were found to 
be compatible with observational data. 

Zhao & Mao (1996) pointed out that except 
for very contrived models, the higher reported op- 
tical depth values are in conflict with the total 
lensable material in the inner Galaxy. The opti- 



cal depth for the low-latitude red clump subsam- 
ple from MACHO Alcock et al. (1996) is between 
(3 - 10) x 10~ 6 at the 95% confidence level. This 
would require contrived models involving the Sun- 
center line being the diagonal line of a homoge- 
neous rectangular bar. This discrepancy has been 
highlighted by Kuijken (1999) and by Binney, Bis- 
santz, & Gerhard (2000) for axisymmetric models. 
Both confirmed that the need of extremely con- 
trived bar and disc models if the optical depth is 
as large as (3 - 10) x 10~ 6 . 

The most recent analysis of bulge microlensing, 
however, suggests a lower optical depth in gen- 
eral. Statistics have improved as we now have 
over 500 events. In particular, the value for red 
clump giants is the range (1.4 — 2) x 10~ 6 depend- 
ing on whether the long duration events are in- 
cluded (Popowski 2001, private communications). 
Sevenster & Kalnajs (2001) suggested that the 3 
kpc arm might correspond to a stellar ring, and 
also account for a large fraction of the microlenses. 

Without making detailed models of the non- 
axisymmetric bar and the ring, we can do a 
quick check of whether our models are consis- 
tent with the microlensing data by scaling the 
results of Zhao, Rich, & Spergel (1996) and Peale 
(1998) for the masses of baryonic components in 
our models. If we use only the p2 bar compo- 
nent (mass 1O 1O M0), then the optical depth is 
- 1O- 6 (M 2 /1O 1O M ) - 10~ 6 (Zhao, Spergel, & 
Rich 1995). In our model A\ there is another 
~ 10 10 M Q of baryons inside 3.5 kpc which is for- 
mally in the ps axisymmetric disk. As we have 
discussed, this thin unperturbed disk inside the 
non-axisymmctric bar is only a mathematical sim- 
plification. It is unlikely that it can stay thin or 
even can survive inside the bar. The most prob- 
able scenario is that the bar was formed from 
the disk at some epoch, when the disk became 
unstable. In this case the disk becomes the bar 
and we should expect that the inner disk has the 
same distribution as the bar. In this case, the mi- 
crolensing counts should include all the baryonic 
mass inside 3.5 kpc, which doubles the estimate 
of the optical depth. Then our best models A\ 
and B\ are consistent with an optical depth of 
(1.4-2) x 10" 6 . 
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5.3. Comparison with other work 

Our estimate for the virial mass of the M31 halo 
1.4 x 10 12 M Q is consistent with earlier estimates 
which used kinematics of satellites. Courteau 
& van den Bergh (1999) found a total mass of 
(1.33 ±0.18) x 10 12 M Q inside 260 kpc, and Evans 
ct al. (2000) and Evans & Wilkinson (2000) found 
(1.23± |) x 10 12 M Q . The main difference between 
our results and Evans et al. (2000) is that they as- 
sumed a halo density profile very different from an 
NFW profile. 

For the Milky Way halo, our models A\ and B\ 
predict a mass A^ fl<10 okpc = (^-8 — 6) x 10 11 M©, 
consistent with what Dehnen & Binncy (1998) 
find in their "standard models 1-4", -^ fl<100 kpc = 
(6.0 — 6.6) x 10 11 M Q . Our disk also has a mass of 
(4 — 5) x 10 10 Mq, consistent with their finding of 
M d = 4.2-5.1 x 10 10 M Q . 

Our results and conclusions differ from those of 
Hernandez, Avila-Reese, & Firmani (2001), who 
model the Milky Way in the framework of a cosmo- 
logical scenario. Hernandez, Avila-Reese, & Fir- 
mani (2001) find that the standard halo profiles 
are not compatible with the observed properties of 
the Milky Way and that "the rotation curve for our 
Galaxy implies the presence of a constant density 
core in its dark matter halo" . The gross features of 
our approach and that of Hernandez, Avila-Reese, 
& Firmani (2001) are similar. The parameters of 
the cosmological models are only slightly different. 
Adiabatic infall corrections are also used, and the 
dark matter profiles are also similar. The main 
difference was in the mass of the halo assumed for 
the Galaxy. The average virial mass in the mod- 
els of Hernandez, Avila-Reese, & Firmani (2001) 
was assumed to be M v ; r = 2.8 x 10 12 M Q , which 
is almost three times larger than the virial mass 
in our favored models A\ and B\. If we ran our 
models for this large virial mass, we would also re- 
ject the model (our model A4 has the largest virial 
mass, M v i r = 2.0 x 10 12 M Q , and even this model 
is marginal, with a rather large rotation velocity 
and an abnormally low concentration). The reason 
that Hernandez, Avila-Reese, & Firmani (2001) 
adopted such a large virial mass is their perhaps 
overly stringent interpretation of the Kochanck 
(1996) estimate of the mass inside 50 kpc. Hernan- 
dez, Avila-Reese, & Firmani (2001) assume that 
the mass was M 50 = (4.9 ± 0.5) x 10 11 M , while 



Kochanek (1996) gives M 50 = (3.2-5.5) xlO n M 
at the 90% confidence level if Leo I is excluded 
from the analysis. For comparison, our models Ai 
and Bi have M 50 = (3.7 - 3.8) x 10 11 M , which 
is compatible with Kochanek (1996). Wilkinson & 
Evans (1999) give even larger uncertainties on the 
mass: M 50 = 5.4±°;| x 1O U M . 

There is no doubt that estimates of the mass 
at large radii, 50 kpc - 100 kpc, are crucial for 
constraining the halo parameters. Values of M50 
larger than 5 x 10 11 M are very likely to be dif- 
ficult to reconcile with the standard cosmological 
models. The problem is that the observational 
constraints are still very much uncertain. 

It is interesting to compare our results with the 
results of Binney & Evans (2001), who came to 
the conclusion that "cuspy haloes favored by the 
Cold Dark Matter cosmology (and its variants) are 
inconsistent with the observational data". Their 
main argument was a combination of three con- 
straints. (1) The mass of baryonic matter in the 
disk and bulge inside the solar radius should be 
larger than 3.9 x 10 10 M to produce enough mi- 
crolensing events. (2) the surface density of the 
dark matter inside 1.1 kpc from the plane at the 
solar radius should be 30 ± 15M pc -2 to satisfy 
the constraints from local steller kinematics and 
stellar and gas content. (3) The sum of all compo- 
nents should reproduce the circular velocity curve 
at 2-4 kpc. Our models A\ and B\ easily sat- 
isfy the first two constraints. For example, there 
is a mass of 4.5 x 10 10 M of baryons inside the 
solar radius in our Model A\, and the dark mat- 
ter surface density is 22M pc~ 2 . In our modeling 
we do not use the circular velocity at 2-4 kpc be- 
cause corrections due to non-circular motions are 
large and uncertain. Still, even if we accept the 
Binney & Evans results for the circular velocity, 
the differences appear to be small — 10 km/s for 
2-4 kpc, which is within the observational uncer- 
tainty. We also note that that the treatment of 
the dark matter and the bar by Binney & Evans 
are not realistic. The profile of the dark matter 
was assumed to be unmodified NFW, which as 
we have discussed is not appropriate for the inner 
part of our Galaxy where adiabatic compression 
will have substantially modified the dark matter 
profile. The bulge was modeled as an elliptical ex- 
ponential disk with scale-length 1 kpc. We have 
used more realistic models for the dark matter and 
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a more detailed bar model which matches numer- 
ous observational constraints. We also rely on the 
more detailed treatment of lensing events made by 
Zhao, Spergel, & Rich (1995) and by Peale (1998) 
for the same bar used in our analysis. These dif- 
ferences explain why we come to different conclu- 
sions. 

6. Conclusions 

We study models of the Milky Way and the An- 
dromeda galaxies based on the standard paradigm 
of disk and dark matter halo formation within a 
ACDM cosmology. The models produced accept- 
able fits for numerous observational data. Here we 
itemize our main conclusions. 

(I) Cuspy NFW profiles may be compatible 
with observational data on MW and M31. 

(II) The compression of the dark matter by 
baryons is a very significant effect in the inner 
parts of galaxies, and must be taken into account. 
It provides a coupling between the baryons and 
the dark matter. The result of this coupling is 
that the total density distribution does not show 
any wiggles or features corresponding to the tran- 
sition from one component to another. Thus, the 
disk-halo "conspiracy" is a natural prediction of 
our models. 

(III) The rate of microlensing events is a strong 
constraining factor for models with standard cuspy 
halo density profiles. These models arc compatible 
with the recent low estimates of the optical depth 
if all the stellar material inside 3.5 kpc is optimally 
distributed in a « 2 x 10 10 M bar. 

(IV) In none of our models does the dark matter 
dominate the central parts (« 3 kpc) of the galax- 
ies. Since both the MW and M31 appear to be 
typical high surface brightness galaxies (HSB), we 
expect that this is true for other HSB galaxies. If 
there has been no transfer of angular momentum 
between the baryons and the dark matter, HSB 
galaxies should have sub-maximal disks. 

(V) Including the effects of a modest transfer of 
angular momentum from baryons to dark matter 
(expected during the non-axisymmetric process of 
bar formation) can produce models with 2-3 times 
less dark matter in the central 3 kpc region of the 
Galaxy. This is likely to allow sustanance of fast 
rotating bars such as those observed in the Milky 
Way and M31. These models also allow a slightly 



heavier disk without producing too high a rotation 
curve. 

(VI) Our dynamical models suggest that the 
virial mass of the dark matter halo hosting the 
Milky Way and M31 is in the range 10 12 M < 
M vir < 2 x 1O 12 M . We use the predicted halo 
mass function to estimate the number density of 
such halos and compare this with the observed 
number density of galaxies. We find that if the K- 
band luminosity of the Milky Way is Mk = — 24, 
as found by Drimmel & Spergel (2001) based on 
direct fits to DIRBE data, then the predicted num- 
ber density of dark matter halos hosting Milky 
Way galaxies is not inconsistent with the observed 
K-band luminosity function of Kochanek et al. 
(2001). If the Milky Way has a larger luminos- 
ity of M K = -24.4 to - 24.7, as would be the 
case if it lies on the Tully Fisher relation found by 
Tully & Pierce (2000) and Rothberg et al. (2000), 
then the number density of "Milky Way" halos is 
more than a factor of two higher than the observed 
number density of galaxies at this luminosity. In- 
creasing the virial mass of the Milky Way's halo 
to the upper limit allowed by the dynamical mod- 
elling (M vir < 2 x 10 12 M Q ) allows the fainter end 
of the TF range to be accomodated, but not the 
brightest values. 

(VII) A significant fraction of the baryons 
within the virial radius of the halo must not be 
in the disk or bulge of the Milky Way and M31. 
All acceptable models required that the disk and 
bulge should contain not more than «l/4-l/2 of 
the baryons expected to be inside the virial radius 
of the halo in the absence of feedback. It is beyond 
the scope of this paper to address the state and the 
location of the "lost" baryons, but this is in keep- 
ing with the usual assumptions of semi-analytic 
models, and we address this issue in Paper II. We 
note that this conclusion is based on two indepen- 
dent arguments. Models which have more than 
half of the baryons in the central luminous part 
fail to produce acceptable fits for either the rota- 
tion curve or for the local surface density. When 
compared with the observed luminosity function, 
these models also produce too large a number den- 
sity of bright galaxies (see (VI) above). This is 
related to the usual "over-cooling" problem (see 
Balogh et al. (2001) for a recent discussion). 
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